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1 Introduction 

Through this paper we consider nearly integrable isochronous Hamiltonian systems as 

2 

H fl =GJ-I + ^- + (cosq-l) + f if(ip,q), (1.1) 

where {<p,q) G T" x T 1 := (R"/27rZ n ) x (R/2ttZ) are the angle variables, (I,p) G R™ x R 1 are the 
action variables and /x > is a small real parameter. Hamiltonian TL^ describes a system of n isochronous 
harmonic oscillators of frequencies u> weakly coupled with a pendulum. 

When fj, = the energy uiili of each oscillator is a constant of the motion. The problem of Arnold's 
diffusion is whether, for /i ^ 0, there exist motions whose net effect is to transfer energy from one oscillator 
to others. This problem has been broadly investigated by many authors also for non-isochronous systems, 
see for example ||, |uj, jnj and p5| . In this paper we focus on isochronous systems for which, in 
order to exclude trivial drifts of the actions due to resonance phenomena, it is standard to assume a 
diophantine condition for the frequency vector u>. Precisely we will always suppose 

• (HI) 3 7 > 0, t > n such that \u ■ k\ > j/\k\ T , Vfc G Z™, k ^ 0. 

The existence of Arnold's diffusion is usually proved following the mechanism proposed in || . First one 
remarks that, for fi = 0, Hamiltonian TL^ admits a continuous family of n-dimensional partially hyperbolic 
invariant tori 7} = {(tp, I. q,p) G T™ x R n x T 1 x R 1 | I = I , q = p = 0} possessing stable and unstable 
manifolds W s (T I(i ) = W u (T Io ) = {((p,I,q,p) G T" x R" x T 1 x R 1 | / = I , p 2 /2 + (cosq - 1) = 0}. 
Then Arnold's mechanism is based on the following three main steps. 

Step (i) To prove that, for n small enough, the perturbed stable and unstable manifolds W^Tj^) and 
Wp(7j*) split and intersect transversally ("splitting of the separatrices"); 

Step (ii) To prove the existence of a chain of tori connected by heteroclinic orbits ("transition chain"); 

Step (Hi) To prove, by a shadowing type argument, the existence of an orbit such that the action 
variables I undergo a variation of 0(1) in a certain time Td called the diffusion time. 

1 Supported by M.U.R.S.T. Variational Methods and Nonlinear Differential Equations. 



We point out that for isochronous systems assumption (HI) implies that all the unperturbed tori 7} , 
with their stable and unstable manifolds, persist, for fi small enough, being just sligthly deformed. For 
this reason the construction of the "transition chain" of step (n) is a straightforward consequence of 
step (i). This also happens, for the peculiar choice of the perturbation, in the non-isochronous system 
considered in On the other hand, this is not the case for general non-isochronous systems where the 
surviving perturbed tori are separated by the gaps appearing in KAM constructions, making the existence 
of chains of tori a difficult matter, see pi]] . We quote paper (2j| for a somewhat different mechanism of 
diffusion where step (n) is bypassed using Mather's theory. 

In the present paper we address, for isochronous systems, the following two main questions 

1) Shadowing theorems and estimates of the diffusion time; 

2) Splitting of separatrices. 

Problem 1) has been intensively studied in the last years, see for example [lCj] , [fTl| , , [Ti|] and pl[ | 
(we underline that JTo| . ]IT[ , |l2|] and 21 deal also with non-isochronous systems). Our general shadowing 



theorem (thm. 2.3-thm.3.i) improves -for isochronous systems- the estimates on the diffusion time obtained 



in the foremcntioncd papers. 

The estimate on the diffusion time that we obtain (see expression (2.23)), once it is verified that the 
stable and the unstable manifolds split, is roughly the following: the diffusion time Td is estimated by 
the product of the number of heteroclinic transitions k (— number of tori forming the transition chain 
= heteroclinic jump/splitting) and of the time T s required for a single transition, namely Td = kT s . 
The time for a single transition T s is bounded by the maximum time between the "ergodization time" 
of the torus T" run by the linear flow cut, and the time needed to "shadow" homoclinic orbits for the 
quasi-periodically forced pendulum. 

In order to highlight the improvement of our estimate of the diffusion time let us consider the par- 
ticular case of "a-priori unstable" systems, i.e. when the frequency vector uj is considered as a constant 
independent of any parameter. In such a case it is easy to evaluate, using the classical Poincare-Melnikov 
theory, that the splitting of the separatrices is O(fi). Then our shadowing theorem yields the estimate 
for the diffusion time Td — 0((l//x) log(l//x)), see thm. 



2.j,-thm.3.3 



Such estimate answers to a question raised in |L9| (sec. 7) proving that, at least for isochronous 
systems, it is possible to reach the maximal speed of diffusion fi/\ log/i|. On the contrary the estimate on 
the diffusion time obtained in Q is Td >> 0(exp (1/^)) and is improved in [[ll] to be Td = 0(exp (l//j)). 
Recently in |l0| by means of Mather's theory the estimate on the diffusion time has been improved to 
be T d = 0(l/^ 2r+1 ). In Jl|l it is obtained via geometric methods that T d = 0(l/jU r+1 ). It is worth 



pointing out that the estimates given in 1 10 and |12| , which yet provide a diffusion time polinomial in 
the splitting, depend on the diophantinc exponent r and hence on the number of rotators n. On the 
contrary our estimate is independent of n. 

The main reason for which we are able to improve also the estimates of |lCj and |L2| is that our 
shadowing orbit can be chosen, at each transition, to approach the homoclinic point, only up to a 
distance 0(1) and not 0(\i) like in jllj and |lS^ |. This implies that the time spent by our diffusion orbit 
at each transition is T s = O (log(l / fi)). Since the number of tori forming the transition chain is equal to 
0(l/splitting) = 0(1/ fi) the diffusion time is finally estimated by Td = 0((l//i) log(l//z)). 

Regarding the method of proof, we use a variational technique inspired by || and [^) . One advantage 
of this approach is that the same arguments can be also used when the hyperbolic part is a general 
Hamiltonian in R 2m , m > 1, possessing one hyperbolic equilibrium and a transversal homoclinic orbit. 
Nevertheless we have developed all the details in the case that the hyperbolic part is the standard 
one-dimensional pendulum because it is the model equation to study Arnold's diffusion near a simple- 
resonance. 



Furthermore we also remark that our proof of theorem 3.2 is completely self-contained in the sense 



that, unlike the known approaches (excepted |25|), we do not make use of any KAM- type result for 



proving, under assumption (HI), the persistence of invariant tori, see thm. 3.1. 

In sections 4-5 we study problem 2). Detecting the splitting of the separatrices becomes a very difficult 



problem when the frequency vector to = u e depends on some small parameter e and contains some "fast 
frequencies" = 0(1 /e b ), b > 0. Indeed, in this case, the oscillations of the Melnikov function along 
some directions turn out to be exponentially small with respect to e and then the naive Poincare-Melnikov 
expansion provides a valid measure of the splitting only for fj, exponentially small with respect to e. Much 
literature in the last years has been devoted to overcome this problem, see for example [f[3| , [p| , |l6| , 
|[20|l and (23). In the present paper, in order to justify the dominance of the Poincare-Melnikov function 
when /i = 0(e p ), we extend the approach originally used in |Q for dealing with rapidly periodic forced 



systems. Up to a change of variables close to the identity, we prove (see thm. ^.1) an exponentially small 



upper bound for the Fourier coefficients of the splitting. As an application we provide some results (thm. 



\5.\ , thm. \5.4 ) on the splitting of the separatrices and the diffusion time (thm. |5.^ ) for three time scales 
systems 

H e = +e a p-h + ^r + (cosq- 1) + fi(cosq - l)f(<p), h e B},I 2 ,P e R n_1 , n > 2, 

V£ 2 

for /ie -3 / 2 small. This improves the main theorem I in [^2| which holds for \i — e p , p > 2 + a. With 
respect to |Tfi| , which deals for more general systems, we remark that our results hold in any dimension, 
while the results of lll| are proved for 2 rotators only. 



Theorem 4.2 is also the starting point for proving the splitting of the separatrices in presence of two 
high frequencies, assuming as in ]i^|,[p6|, p~of and |||] suitable hypotheses on the perturbation term. We 
do not address this problem in this paper. 



On the other hand theorem 5.1 is the starting point to prove, for n > 3, the existence of diffusion 
solutions such that the action variables I2 undergo a variation O(l) in polynomial time (while the I\ 
action variable does not change considerably). This phenomenon can not be deduced by the estimates, 
given in jlfj and |22| , on the "determinant of the splitting" which does not distinguish among "slow" 
and "fast" directions and would give rise to exponentially large diffusion times. This type of results are 
contained in the forthcoming paper || . 

The paper is organized as follows: in section 2 we prove the shadowing theorem when the perturbation 
term is /(</?, q) = (1 — cos q)f(ip). In section 3 we show how to prove the theorem for general perturbation 
terms f(ip,q). In section 4 we provide the theorem on the Fourier coefficients of the splitting and in 
section 5 we consider three time scales systems. 

The results of this paper have been announced in Q . 

After this paper was completed we learned by prof. Bolotin about the recent preprint which deals 
with a-priori unstable Hamiltonian systems time periodically forced. Among many results, in theorem 
2 of p4| a shadowing theorem for a symplectic separatrix map is proved providing an estimate on the 
diffusion speed log/z|, like our. However theorem 2 applies just to an approximation of the symplectic 
map describing the true dynamics of the system and, moreover, requires the hyperbolic part to be just 
two dimensional. 

Acknowledgments: The first author wishes to thank Prof. G. Gallavotti for stimulating discussions. 



2 The shadowing theorem 

We first develop our approach when the perturbation term f(ip, q) = (1 — cosq)f(ip) so that the tori 7} 
are still invariant for fi ^ 0. The equations of motion derived by Hamiltonian TL^ are 

ip = w, I = -//(l - cosq) d v f(ip), q = P, p = sinq - [isinq f((p). (2.1) 

The dynamics on the angles tp is given by <p(t) = cut + A so that (|]l]) are reduced to the quasi-periodically 
forced pendulum 

— q + sing = fi sin 5 f(ojt + A) (2-2) 

corresponding to the Lagrangian 

C^q, q,t) = ^- + (l- cosq) + ^(cosq - l)f(ut + A). (2.3) 



For each solution q(t) of (2.2) one recovers the dynamics of the actions I(t) by quadratures in (2.1). 



2.1 1-bump homoclinic and heteroclinic solutions 



For /x = equation (2.2) is autonomous and possesses the homoclinic (mod. 2-7r) solutions qe(t) = 
4 arctg(exp (t — 9)), 9 e R. Using the Contraction Mapping Theorem we now prove that, near the 
unperturbed homoc linic solutions qe(t), there exist, fo r /i small enough, "pseudo-homoclinic solutions" 
1a e W °f equation (2^). q A g (t) are true solutions of (2^) in (— oo, 9) and (9, +oo); at time t = 9 such 
pseudo-solutions are glued with continuity at value q A g (9) = n and for t — > ±oo are asymptotic to the 
equilibrium mod 2ir. 

Lemma 2.1 There exist Ho,Cq > such that VO < /i < (1q, Vw G R", V# G R, i/iere exists a unique 
function q A g (t) : R — > R, smooth in (A, 9, /i), such that 



i) q A g(t) is a solution of (2.i ) in each interval (— oo, 9) and (9, +oo) and q A Jff) = n; 
it) max [\q^ j0 (t) - g e (t)|, \q^ e (t) - < Co/xexp(--^) ; Vi G R; 

Hi) q%(t)=q^ +k2Kje (t), VfcGZ"; 

iv ) llfi+nit + V) = Qa+wvA^' V61 ' ^ e R - 

v) max (|0 A <£ ffl (t)|, \d A q% e (t)\, \u ■ d A q%{t)\, \u ■ d A q%{t)\) = o(/xexp(-J^)) . 

Proof. Proof in the appendix. Note that in (v) the bound of \u> ■ d A q A gitylA^ ' ^4<?a e(*)l ^ s uniform 



We can then define the function : T" x R — > R as the action functional of Lagrangian (2.3) 
evaluated on the "1-bump pseudo-homoclinic solutions" q A g(t), namely 



+00 



^ (<?%(*), dt+ / cM g (t),fig(t),t) dt 



and the "homoclinic function" G M : T™ — > R as 

G I1 {A)=F I1 {A,Q). 



(2.4) 



(2.5) 



Since q A g (t) converges exponentially fast to 0, mod 2ir, the integrals in (2.4) are convergent. Note that 
the homoclinic function G M is independent of Iq. By property (v) of lemma 2.1 the following invariancc 
property holds 

F^A,9 + n) = F^A + Lun,9), W, rj £ R, 

and in particular 

F^(A, 6) = Gp(A + ui6), V0 g R. (2.6) 

Remark 2.1 The homoclinic function G M is the difference between the generating functions Io (A, qo) 
of the stable and the unstable manifolds W*' w (Tj ) (which in this case are exact Lagrangian manifolds) 
at section qo — n, namely G^(A) — S~ I(j (A, ir) — S~T j (A, it) . Indeed can be easily verified that 



r+oo 

S+ Jo (A,gb)=VA- / 
Jo 



+ (1 - cos + n(cos(ft„(t) - !)/(«* + A) dt, 



where q Aqo (t) is the unique solution of (2.i) near qo(t) with q A qo (0) = qo and lim t ^ +00 q Aqg (t) = 2ir. 
Analogously 



S „,i ( A ><}o) :=Io-A + 



(&*,(*)) 



+ (1 - cos<# (t)) + /x(cos< ;9o (t) - l)/(wt + A) dt 



where q A qo (t) is the unique solution of (2.i) near qo(t) with q A qo (0) — qo and lim^-oo q A _ (i) = 0. 



Lemma 2.2 The derivative of 9 — > F fl (A,9) satisfies 



Proof. There holds 



deF M (A,6) = 



+ 00 



(2.7) 



i%(t)deq%{t) + (sinq%(t) - [j, sin q^ e (t)f(wt + A)^d e q^ e (t) dt. 
Integrating by parts and using also that q\ g (t) solves (|2.2| ), we obtain 



dgF^A,9) = \{q%o) 2 {0-) \{i%o) 2 (0 + ) + \deq%{t)q% e {t) 



d eQA.e(t)qA,9(t) 



(2.8) 



Since \/9 S R q\ g (8) = 7r, deriving in we obtain deq A g (9) + q A g (9) = 0; hence from (2.8) and using 
that lim^-i-oo q A g (t) = we deduce lemma 2.2. ■ 

By lemma |2.2| if dgF fJ/ (A, 9) = then q A g (t) is a true homoclinic (mod. 2tt) solution of (2.2). Then, 
for each I £ R", 

(wt + AJ^t),q%{t),q%{t) 



where 



In{t) =Jo — fi (1 - cosg^ e (s))d v f(ujs + A)ds 



(2.9) 



(2.10) 



is a solution of 7i M emanating at t = — 00 from torus 7} . Since q A g converges exponentially fast to the 
equilibrium, the "jump" in the action variables I M (+oo) — 1$ is finite. We shall speak of homoclinic orbit 
to the torus 7} when the jump is zero, and of heteroclinic from 7} to Tj- (+00) when the jump is not 
zero. Moreover the next lemma says that such jump is given by d^F^A, 9): 



Lemma 2.3 Let dgF fl (A,9) = then I^if) given in ( 2.1C ) satisfies 

/*+oo 



d A FJA, 



I^t) dt = J M (+oo) - I < +00. 



(2.11) 



In particular if (A, 9) is a critical point of F^{A,9) then (2.S) in a homoclinic orbit to torus Tj . 
Proof. There holds 

/+00 
iA fi (t)dAq%(t) + smq% e (t)d A q%(t) 
-OO 

- ^sinq^ e (t)f(ujt + A)d A qA, e ( t ) - M 1 - C0S 9A,e(*))<WM + A ) dt - 
Integrating by parts, since q A g (t) solves (p^), and using that lim^±oo q A g (t) = 0, we deduce 

r+oc 



/+00 
-,i(l - cos q^ g (t))d v f(ujt + A) dt. 
-OO 



(2.12) 



We deduce from ( 2.12 ) equality ( 2.11 ) 
By the invariance property ( 



2.6) if B is a critical point of the homoclinic function G^, then, for all 



(A, 9) such that A + uj9 = B, (2.9) are homoclinic solutions to each torus Tj . These homoclinics are 
not geometrically distinct since, by the autonomy of Ti^, they are all obtained by time translation of the 
same homoclinic orbit. By the Lusternik-Schirelman category theory, since cat T" = n + 1, the function 
Gfj, : T" — > R has at least n + 1 distinct critical points. This proves (see also [pof ) 



Theorem 2.1 Let < p, < po. V/o S R" there exist at least n + 1 homoclinic orbits geometrically 
dinstict to Ti a . 



From the conservation of energy a heteroclinic orbit between 7} and 7}^ , if any, must satisfy the 
energy relation 

uj ■ I = co ■ I' . (2.13) 



By lemma |2J5| a critical point of F^ IoJ ^(A,9), defined by F^j^A, 6) = F^(A,9) - (I' — I ) ■ A = 
G /J ,(A + lu O) — ( Iq — I ) ■ A, gives rise to a heteroclinic solution joining the tori 7/ to 7}^. If the energy 
condition (2.13) holds then the function F^ ^ j> (A, 9) satisfies the invariance property 

F^i o j,(A,d)=G^A + oj0)-(l{ ) -I Q )-(A + uj9)=G, l j o j,(A + u8). (2.14) 

where 

G M ,{B) := G^(B) - (I' - I ) ■ B. (2.15) 

Note that G^jgj^ is not 27rZ"-periodic, and it might possess no critical point even for |/q — Jo I small. 
However near a homoclinic orbit to 7} satisfying some "transversality condition" there exist heteroclinic 
solutions connecting nearby tori 7}' . As an example, the following theorem holds, where B p (A ) denotes 



an open ball in R" (covering space of T"). 

Theorem 2.2 Assume that there exist Aq € T", S > and p > such that inigs (A ) G> > m fs p (A ) G fJ- + 
5. Then for all IqtI'o £ R-™ satisfying (Iq — Iq) ■ to — and \Iq — I Q \ < S/(2p) there exists a heteroclinic 
solution of Tl^ connecting Ti toTj^. 

2.2 The fc-bump pseudo-homoclinic solutions 

We prove in the next lemma the existence of pseudo-homoclinic solutions q\ g (t) of the quasi-periodically 
forced pendulum {2/2) which turn k times along the separatrices and are asymptotic to the equilibrium 
for t — > ±oo. Such pseudo-homoclinics q\ g (t) are found, via the Contraction Mapping The orem, as small 



perturbations of a chain of "1-bump pseudo-homoclinic solutions" obtained in lemma 2.1 



Lemma 2.4 There exist C 1: L x > such that Vuj G R", VO < /i < p , Vfc e N, VL > L x , V6»i < . . . < 9 k 
with mini — $i| > L, there exists a unique pseudo-homoclinic solution q\ g (t) : R — > R, smooth in 
{A, 9, p) which is a true solution of (2.i ) in each interval (— oo, 9\), (9i, 9i + \) (i — 1, . . . , k — 1), (9 k , +oo) 
and 

=7r(2t-l), q^ d (t) = q% 1 (t) in (-00,61) andq^ e (t) = 2n{k-l)+q% k {t) in (9 k ,+oo); 

") lki,e-9A9 < lk I -(Ji) < Ciexp(-CiL) where Ji = (6i, (9 t + 9 l+ i)/2), V i = 1, . . . , k - 1; 



Q'A.e-^A+iH^ 1 ' 00 ^) - Ci ex P(-Ci^) where J[ = ((0* + <+ i)/2, i+ i), V i = 1, . . . , k- 1; 



«A,fl-H/(* + »») = «A+wij,fl(*)» V0,77eR. 



Proof. In the appendix. ■ 

We co nside r the Lagrangian action functional evaluated on the pseudo-homoclinic solutions q\ g given 
by lemma 2A depending on n + k variables 

F^A 1 ,...,A n ,8 1 ,...,9 k ) = [ 1 C^(q L Afi (t),q L Afi (t),t) dt + 



3C 



rOi+i f-\-oo 

E/ ^(<lAAt),qA,e(t),t)dt+ C^q L A£ {t) 1 q L A ,e{t),t) dt 



By lemma 2.4 -v the following invariance property holds 

F* (A, 9 + 7!)= F* (A + r,u, 9) , V0, n e R. 
Let : T" x R fc — * R be the "fc-bump heteroclinic function" defined by 

J*{A,tf) :=F*{A,0)-{I' Q -I Q )-A. 

Arguing as in lemma |2.3| we have 



(2.16) 
(2.17) 

Lemma 2.5 Vib,io <= R "> */ (A 0) is a critical point of T*{A,0), then (cut + A, I^{t), q\ d (t), q\ e (t)) 
where I^(t) = Iq — A* (1 — cos 5^ g(s))d v> f(u>s + A)ds is a heteroclinic solution connecting 7} n to . 



By lemma 2J5 we need to find critical points of F^{A, 8). When mini(^ + i — 9A — ► +00 the "k- 
bump homoclinic function" F*(A,6) turns out to be well approximated simply by the sum of F^{A,9i) 
according to the following lemma. We set 6*o = — 00 and 9^+i = +00. 

Lemma 2.6 There exist C 2 ,L 2 > such that G R", VO < (X < /i , ML > L 2 , V6>i < . . . < 9 k with 
minj(0j + i - OA > L 



with 



F*(A, 1 ,...,0 k )=J2 F^A, 9,) + A Bi-iA, i+ i) 

i=l i=l 

\Riifx, A, O^A, < C2 exp(-C 2 L). 



(2.18) 



Proof. We can write 
F*{A,0 x ,...,0 k ) = 



(e 1 +e 2 )/2 



fc-i 

E 

i=2 



(Si-i+9i)/2 



9i+B i+1 )/2 



) fc _!+e fc )/2 



r+00 

^(3A,fl(*)»3A,fl(*)»<)+ / A.(3A,fl(*)>9A,9(*)'*)J- 



We define 

Rj~(n,A,Oi-i,0i) 

R+(n,A,0 u i+1 ) - 



(e,+e I+1 )/2 



^((«,fl(*).«A,fl(t),*)dt- / £^A A Xt),q% (t),t) dt, 



+ OO 



^(3A,fl(*)»«,fl(*)»*)*- / A.(<fl 4 (*).9W*).*)d* 



where e . is the 1-bump pseudo-homoclinic solution obtained in lemma 2.1. Recalling the definition 2.4 
of F^A, 9) we have 

F^(A,8 1 ,...,0 k ) = F^A,9 1 ) + R+(n,A,9 1 ,9 2 ) 
fc-i 



^ F M (A, 0;) + A, 0;_ l5 0;) + R+(fx, A, 0i,0 i+ i) 



i=2 



+ F^A^kj+Rl^AQk-iJk)- 



Setting Ri — R^ + Rf we derive the expression (2.18). In order to complete the proof, it is enough to 
show the existence of C2, L 2 > such that Vw £ R", for all < fx < fio, VX > L 2 , V#i < . . . < 0u with 
min, (^ + i — 9{) > L, for all i = 1, . . . , k 



\Rf(fx,A,0i,0 i+1 )\ <C 2 exp(-C 2 L). 



(2.19) 



We write the proof for Rj . We have 

•{9i+e i+x )/2 



J9i 



+00 

h+e i+1 )/2 



(£ M (q% e (t) ,qi (t),t)- £ M {(frfit (*) (*) . *) * 



(2.20) 



By lemma ^^-(m) the homoclinic orbit satisfies max(|g^ e .(i)|, e .(t)|) < Cexp(— |t — 0j|/2). Hence, 
for all #i < . . . < 9k with min^^+i — 9i) > L, 



+ OC 



i+e i+ i)/2 

From lemma [Q|-(m) we also deduce that 

0i+e i+1 )/2 



(2.21) 



(2.22) 



From (|2.20|) , (p.2l|) and (|2.22j) we deduce fl2.19| ) and hence the lemma 



2.3 The diffusion orbit 

We are now able to consider the existence of the shadowing orbit. We give an example of condition on 
which implies the existence of diffusion orbits. 

Condition 2.1 ("Splitting condition") There exist Aq € T™, <5 > 0, < a < p such that 

• (i) ™f dBp{Ao) Gf, > inf Bp(Ao) G M + 5; 

• (H) sup Ba(Ao) G p < I + inf Bp(Ao) G^; 

. (m) d({A e B P (A ) I Gfj, (A) < 5/2+inf Sp(Ao) G»},{A e B P (A ) \ G M (A) > 36/4+mf Bp(Ao) G p }) > 
2a. 

Remark 2.2 IfG^ possesses a non-degenerate minimum in Aq the "splitting condition" above is satisfied, 
for p sufficiently small, choosing 5 = (minAj)/3 2 /4 and a = (p/8)^/(mhxi Ai)/(maxi Ai) where Ai are the 
positive eigenvalues of D 2 G^(Aq). 



Remark 2.3 _B p (j4o), open ball of radius p in R™ (the covering space ofT n ), could be replaced by a 
bounded open subset U of R™ . 

The following shadowing type theorem holds 



Theorem 2.3 Assume (HI) and the "splitting condition" 2.1 . Then \/Io,Iq with uj-Iq — uj-Iq, there is a 
heteroclinic orbit connecting the invariant tori 7} Q and Tp . Moreover there exists C3 > such that V77 > 
small enough the "diffusion time" Td needed to go from a rj-neighbourhood of '7} to a 77 -neighbourhood of 
Tj' g is bounded by 

Hi-iMp max (| hl S \ , J-) + C 3 1 Info) I • 



T d < C 3 



(2.23) 



Remark 2.4 The meaning of (2.2c ) is the following: the diffusion time Td is estimated by the product of 
the number of heteroclinic transitions k — ( heteroclinic jump / splitting ) — — Iq\/6, and of the time 
T s required for a single transition, that is Td = k ■ T s . The time for a single transition T s is bounded by 
the maximum time between the "ergodization time" (l/~fa T ), i.e. the time needed for the flow tot to make 
an a-net of the torus, and the time | In (5 1 needed to "shadow" homoclinic orbits for the forced pendulum 
equation. We use here that these homoclinic orbits are exponentially asymptotic to the equilibrium. 



Remark 2.5 The following proof works if G u p ossesses a local maximum which satisfies a non-degeneracy 
type condition like the "splitting condition" 2.1, while in the approaches developed in fcldll and J§<^/ ; based 
on Mather's theory, diffusion orbits are always built from local minima ofG^. The proof of the shadowing 
theorem when the homoclimic point Aq is a saddle point requires slightly different arguments. For example 
it holds assuming as in J7^/ the condition D 2 G fl (Ao)ui ■ uj ^ 0. 

Proof. Assume with no loss of generality that Aq — and mf Bp ( ) G^{A) — 0. Let us choose the 
number of bumps k as 

-2A-p-\I' Q -I \ 



1. 



(2.24) 



By lemma |2.4| -(?) and lemma 2A-(ii), the trajectory converges exponentially fast to 7} (resp. 7}<) as 
t — » — oo (resp. +oo) from 9\ (resp. 9k). Therefore it is enough to pro ve th e existence of a critical point 
(A, 9) e T n x R fc of the /c-bump heteroclinic function J 7 * defined in ( 2.17 ), such that for some positive 
constant K\ 

1 



_0J <K ^ Iq f J °' pmax(|ln<5|,— 



More precisely we shall enforce 



K 2 \\n5\ < \9 i+1 -9i\ < K 3 max.(\hi6\, 



1 

7a 7 



Vi = 1, 



(2.25) 



(2.26) 



for some positive constants K2,K%. Let (Oi, • • • , f2 n ) be an orthonormal basis of R™ where 

fii = -r— and VL 2 = 7- 

M \r -Io\ 



We recall that u ■ (Iq — Jo) = 0. In order to find a critical point of J 7 ^ we introduce suitable coordinates 
(a 1 ,...,a n ,si,...,Sk) G R" x (~p,p) k defined by 



Vi = 1, . 



where r\i are constants to be chosen later. In these new coordinates the heteroclinic function defined in 

rji + «i - ai n k + s k - a 1 ' 



(2.17) is given by 



^(01,02, ■ ■ ■ ,a„,si, . . . ,s k ) = Fn^Tajtoj, ■ 



-%-I \a 2 . (2.27) 



Using the invariance property ( 2.16|) we see that T k does not depend on the new variable ai : 



•^(ai, «2, ■ • ■ ,a„, si, . . . , Sk) 



F, 



(E 

i=2 



Ojflj, 



K-I0W2 



^(0,a 2 , . . . ,a„,si, . . .,Sfc). 



For simplicity of notation we will still denote Tz{a,2, . . . , a„, s%, . . . , s&) := J-^(0, 02, ... , a n , si, . . . , Sfe). 
We now choose the contants r)i- Let 



M 
C 2 



In 



+ 2/9, 



(2.28) 



where C2 is the constant appearing in lemma 2.6. We shall use the following fact (see |4|): there is C > 
such that, for all intervals J C R of length greater or equal to C/(^a T ), there is 9 € J such that 



d{6u,2iiZ n ) < a. 



(2.29) 



By Q2.29D there is (771, . . . , 77*.) G R fe such that 

n 

77,^1 = Xi, mod27rZ n , |x»| < « and Xi ' ^1 = 0, i.e. Xi = Z Xi ' (2.30) 



3=2 



/ CM 



(2.31) 



By (2.28), (2.31), since S (—p,p) we have that — 0j > ln( 2 ^ a )|; hence, by lemma 2.6, setting 

i?i = Ri(a,2, ■ ■ ■ , a n , Si-i, Si, s;+i) 

R /\r^ g»— 1 + V1-1 - ai Si + rji - ai s l+1 + rjj+i - ai 
- '''IZ. | W) > | W | > | W | 



we get 



|-Ri(a 2 , . . . , a„, Sj_i, s i+ i)\ < — . 



(2.32) 



By lemma 2.6, the invariance property ( |2.16 ), (2.3C) and since G M is 2-7rZ n -periodic, we have 



J 7 ^(a 2 ,...,a n ,si,...,s k ) = ^F^(^^ajQj,^r^j +Ri-\Io- I \a 2 



i=l j=2 
k 



i=l j=2 
n 

i=l j=2 
fe 

Z G ^( a2 + X'^ 2 ' ■ ■ ■ ,a n + Xi,n, Si) + Ri - \I' - h\a2 



o|a 2 



i=i 

where G M (a2, ■ ■ ■ , a n , s) = G M (^" =2 ajilj+sQi). Since the basis (Qi, . . . , f2 n ) is orthonormal the function 
G M satisfies the same properties as G^, i.e. 

sup G M < (J/4, inf G M > S and d({x <= B p (Q) | G M (x) < 5/2}, {x € BJO) \ G^x) > 35/4}) > 2a. 

B«(0) d,B p (0) 



We shall find a critical point of J 7 ^ in 



W 



|(a 2 , . . . ,a n ,s) 



€ R™^ 1 x R* 



(02, • • • , a„, Si) S Sp(0), Vi = 1, . . . , k j. 



TV\ attains its minimum over W at some point (a, s). Notice that by (2.32) 

k 



inf 5* < £*(0, 0) = V G M (0, x.) + fc 



24' 



Since < a for all i = 1, . . . , fe and sup Bct ( ) G> — 'V'^ we nave 

inf ^< 4+4=4. 

4 24 24 



(2.33) 



The theorem is proved if we show that (a, s) € W . Arguing by contradiction assume that (a, s) € 9W. 
Then there is some I € {1, • • • , fc} such that (a + s/) G dB p (0), so that G M (a + x/, s;) > (5. We now 
prove that (a + *{,*);* € (— p, p)} D B p (0) C Z := {x e B p (0) \ G p (x) > 3<5/4}. Indeed, if not, by (|3|), 
for some t £ (— p, p) such that (a + Xht) & B p (0), 



Fu(a,si, ■ ■ ■ ,si-i,t,si + i 



,s k ) < rf(a,s) + (G p (a + X i,t)-G p {a + xi,si))) 

+ \Ri-i{a,si-2,si-i,si) - Ri-i(a,si-2,si-i,t)\ 

+ \Rl(a, S/_i,t,S/) - i? ; (a, S;_i, S;,S /+ i)| 

+ \Ri+i(a,si,si + i,t) - Ri + i(a,si,si + i,Si +2 )\ 



66 
24 



which is wrong since (a, s) is the minimum of jRJ over W. We deduce in particular that, for all i, 

(fi + Xi-,~i) £ ^ U -Sp(0) c - Now, as G M > 3(5/4 in a neighbourhood of dB p (0), our splitting condition 
implies that 

d{{x e B p (0) | G p (x) > 3(5/4} U B p (0) c , {x € B„(0) | G M (s) < <5/2}) > 2a. 
We derive by (|3|) that for all i, 



G^a + XhSi) > 5/2. 



(2.34) 



As a consequence, noting that, from (2.24), |/q — Iq\ p < (fc(5)/24, we deduce that 

2k r- — \ \ j 5 , (5 ,5 ,10(5 7(5 

J 7 ,, (a, s) > fc k k — = k > k — , 

^ y ' ' ~ 2 24 24 24 24' 

contradicting Q2.33 ). The proof of the theorem is complete. ■ 

When the frequency vector w is considered as a constant, independent of any parameter ("a priori- 
unstable case") it is easy to justify the splitting condition 2.1 using the first-order approximation given 
by the Poincare-Melnikov primitive. With a Taylor expansion in p we can easily prove that for p small 
enough 

G P {B) = Const + pT(B) + 0(p 2 ), VB € T n , 
where T : T™ — > R is nothing but the Poincare-Melnikov primitive 



T(B) 



(1 - cos q (t))f (cot + B) dt. 



R 



Hence, if T possesses a proper minimum (resp. maximum) in Aq g R™, i.e 3r > such that inf gB^tAn) -T > 
T(A ) (resp. s\xp SBr ( Ao \T < T(A )) then, for p small enough, the "splitting" condition 2A holds with 
<5 = O(p), p = 0(1) and a — O(l). We re mar k that the previous B r (Ao) could be replaced by a bounded 
open subset U of R". Applying theorem 2.3 we deduce 

Theorem 2.4 Assume (HI) and let T possess a proper minimum (or maximum) Aq, i.e. 3r > such 
that infas r (A )r > T(Aq). Then, for p small enough, the same statement of theorem 2J. holds with a 
diffusion time Td = 0((l/p) \og(\/p)). 



3 More general perturbation terms 

In this section we show how to adapt the arguments of the previous section when dealing with a more gen- 
eral perturbation term f(<p, q). Regarding regularity it is sufficient to have finite large enough smoothness 
for /. The equation of motion derived by Hamiltonian TL^ are 

(p = uj, I = -fj,d 9 f(ip,q), q = p, p = sinq- fid q f(ip,q), (3.1) 
corresponding to the quasi-periodically forced pendulum 

- q + sinq = fi d q f(ujt + A,q). (3.2) 

3.1 Invariant tori in the perturbed system 

The first step is to prove the persistence of invariant tori for /i ^ small enough. It appears that no 
more than the standard Implicit Function Theorem is required to prove the following well known result 
(see for example Jl|] for a proof) 

Theorem 3.1 Letuj satisfy (HI). For h small enough and \/ la G R" systemH^ possesses n-dimensional 
invariant tori Tj q « 7} of the form 

T£ = {/ = Jo + a"(V0, <P = 1>, Q = Q"W0, P = P"W>), V> G T"}, (3.3) 

with Q M (-), -P^(') = 0([i), a^(ip) = 0(h)- Moreover the dynamics on Tj is conjugated to the rotation of 
u for ip. 



We first determine the functions Q^(-), P^() in (3.3). Using the standard Implicit Function Theo- 
rem we prove that there exists a unique quasi-periodic solution q^(t) for the quasi-periodically forced 



pendulum (3.2) which bifurcates from the hyperbolic equilibrium 0. 
Lemma 3.1 Let f € C'(T" x T). For fi small enough there exists a unique quasi-periodic solution 



(q^(t),p^(t)) of (3.i ) with (q^(t) , p^(t)) — 0(h), C 1 -smooth in A. More precisely there exist functions 



Qn^pn : T" ^ R of class C 1 ' 1 , such that (q^(t) , p£(t)) = (Q^(ujt + A), P^(ujt + A)). 

Proof. Let L be the Green operator of the differential operator h — > —D 2 h + h with Dirichlet boundary 
conditions at ±oo. L is explicitely given by L(f) = J R e~' t ~ s 'f(s) ds/2. It results that L is a continuous 
linear operator in the Banach space of the continuous bounded functions from R to R, which we shall 
denote by E. We consider the non-linear operator S:Rx T" x E — > E 

S(h, A, q) :=q- L(q - sing) - [iL(d q f(wt + A, q)). 

S is of class C' -1 . We are looking for a solution q^ of S(n,A,q) = 0. Since S(0,A,0) = and 
d q S(0, A, 0) = Id, by the Implicit Function Theorem there exists, for fi small enough, a unique solution 



q^ A = 0(h)- By (3.2) q A £ C /+1 (R); moreover it is C l 1 -smooth in A. We define the C l 1 -maps 
Q^(-),-P M (-) : T n -> R by 

Q»(A) := g£(0), P»(A) := q%(0). 
By uniqueness we deduce that q^(s + 1) = 9^ +ws (t), Vs, t G R. For t = this yields 

#(«)=<#W0):=Q''(A + W «) and P*X( S ) = Pa+us(®) := P^(A + ws), Vs e R 
proving the lemma.B 



We now define the functions a»(i/)) of (3.3). We impose that (ujt + A,I + a^(ujt + A),Q»(ut + 



A), pv(ojt+A)) satisfy the equations of motions (|3.l|); hence the functions a M (V>) must satisfy the following 
system of equations 

(wVKW^s'm where tf ty) := -(V^/)(V, Q M (V))- (3-4) 



In order to solve (3.4) we expand in Fourier series the functions a AI (-0) = X]fcez™ a k^ lk , 
S/cez™ 9ke lk '^ ■ Each Fourier coefficient Ofe must then satisfy 



i[k ■ uj)ak — jJ.gk, Vfc £ Z r 



(3.5) 



It is necessary for the existence of a solution that g = L, n g ,i (ip)d'ip — 0. This property can be checked 
directly, that is 



Lemma 3.2 We have 



Proof. For all i = 1, 



(3.6) 



(3.7) 



Since (q^(t) , p^(t)) satisfies the pendulum equation ^^l^j^jQ^W = ^(VOj ^'^-^(VO 
sinQ^ip) — nd q f(ip, Q^(ip)) and we deduce that 



We now prove that 



d (P^)) 2 
dV^ 2 



5^Wi(fl0 < (Q"^P' t ) - fy, (Q^P") 



Indeed 



(3.8) 



(3.9) 



3 = 1 



d (P^(V>)) 2 



d-^i 



3¥* 



From (|3.7|), (|3.8| ), fl3.9| ) we finally obtain that 

d^MQ") = ^(/MO") + -cosQ^ + - —J + -^^-^(Q^P^) -^(Q^P")) 



from which property (3.6) follows. ■ 

Since u> satisfies (HI) and / is sufficiently smooth the function defined by 



feeZ",fc^o 



(3.10) 



which formally solves equation (3.5), is well defined and smooth. Indeed since / <G C l the function g* 1 
defined in Q) is C'" 1 and there exists M > such that \g k \ < Af/|fc| i_1 , Vfc G Z™, k ^ 0. By (ill) it 
follows that |a fc | < M/\k\'- x \w • fc| < M|fc| T /(7|fc|'- 1 ). The proof of theorem O is complete. 



3.2 The new symplectic coordinates 

In order to reduce to the previous case we want to put the tori Tj 1 at the origin by a symplectic change 
of variables. Recalling that the tori Tf are isotropic submanifolds we can prove the following lemma 

Lemma 3.3 The transformation of coordinates (J,ip,u,v) — > (I,ip,q,p) defined on the covering space 

R 2(«+l) Q j T n x R" X T x R &y 

I = a"(ip)+ud i) Pi"(ip)-vd, ll Q l "(ip) + J, (p = ip, q = Q t *{iP) + u, p = P f *(i/j)+v (3.11) 
is symplectic. 

PROOF. Set dl A dip = Y%=i dh A d(fi and dJ Adip = Y%=i dJi A dipi. We have 

n 

dlAdip + dpAdq = da^(ip) A dipi + d^d^P^ {ip)") A dipi — dfvd^Q^^ip)) A dtpi 

i=l 

+ dJ AdiP + dv Adu + dP^iijj) A dQ»{ip) + dP»(ip) A du + dv A dQ»(ip). 

Using that the tori Tj 1 are isotropic, that is Y^h=i A dipi + dP^(ip) A dQ^(ip) = 0, and noticing 

that J2 itj ^^/"(Wj A dipi = = Y, itj vE^^Q^dipj A dipi we deduce 

n 

dlAdif + dpAdq = ^ d^d^P^ 1 (ip)) A dipi - divd^Q^ (ip)) A dipi 
i=i 

+ dJ A dip + dv A du + dP^(ip) Adu + dv A dQ^iip) 

n 

= dJ A dip + dv A du + y~] d^P^(ip)du A dipi - d^Q^(ip)dv A dipi 

4=1 

+ dP» (ip) Adu + dv A dQ^ (tp) 
= dJ A dip + dv A du, 



and the transformation (3.11) is symplectic. ■ 

In the new coordinates each invariant torus Tj^ is simply described by {J = Iq, ip S T™, u = v = 0} 
and the new Hamiltonian writes 

v 2 

(/C M ) JCn = E,, + u ■ J + — + (cos u - 1) + Po(M> i>) 

where 

P (/x, tt,^) = ( cos(Q^ + u) - cosQ^ 1 + (sin Q»)u + 1 - cosu) + fi(f(ip, Q^ + u)- f(iP, Q^) - 8 q f(ip, Q*)u 

and Ef,, is the energy of the perturbed invariant torus Tq = {(a^(ip) : ip, Q^(ip), P^iip)); ip £ T"}. Hamil- 
tonian (JCfi) corresponds to the quasi-periodically forced pendulum equation 

— ii + sinu = d u Po(fi,u,ujt + A). (3-12) 

of Lagrangian 

u 2 

= — + (1 - cosu) - P (n,u,wt + A). (3.13) 
Since the Hamiltonian /C M is no more periodic in the variable u we can not directly apply theorem 



we 



2.3 and the arguments of the previous sections require some modifications. Arguing as in lemma 2.1 
deduce that, there exists, for /i small enough, a unique 1-bump pseudo-homoclinic solution g(t), true 



solution of (3.12) in (— oo, 9), (9, +oo), satisfying all the properties of lemma 2.1. Then we define the 
function T a ■ T" x R -> R as 



2 



+ (1 - cosu^ g ) — Po(n, u A g , tut + A) dt 
+ (1 - cos u'X e ) - Pi (p, u A e , wt + A) dt + 2irq% (9) , 



where, Vi € Z, we have set 

Pi(p, u, tot + A) = (^cos(q A (t) + u) — cosq A (t) + sinq A (t) (u — 2ni) + 1 — cosu) 

+ / i(/(w* + A,^(t)+ti)-/(w* + A > ^(i))-(a,/)M + A > ^(i)) («-27ri)). 

Since e converges exponentially fast to for t — > — oo and to 27r for t — > +oo the above integrals are 
convergent. The term 2nq A (9) takes into account that the stable and the unstable manifolds of the tori 
Tj^ are not exact Lagrangian manifolds, see p0| . We define the "homoclinic function" : T™ — > R as 



(3.14) 



It holds also T^A, 9) = Qp,(A + u>9), V# 6 R. Arguing as in lemma [2.4| we can prove the existence of 
fc-bump pseudo-homoclinic solutions u AS , which is a true solution of ( 3.12 ) in eac h interval (— oo,#i), 
(9i, 0i+i) (i = 1, . • • , k — 1), (9k, +00), and satisfying all the properties of lemma 2.4. Then we define the 
"fc-bump heteroclinic function" 



F k JA,6 x ,...,6 k ) 



+ (1 - cosu^ g ) - P (iJL,u\ e ,wt + A) dt + 2wq^(9 1 ) 



' — OO 
fc-1 



E 



+1 (u L A,e? 
2 

+00 („-.L \2 



+ (1 - cosui e ) - Pi{fx, u h Afl , ojt + A) dt + 2irfo(9 i+1 ) 



+ (1 - costing) - P k (iJ,, u£ e ,ut + A) dt- (I' -I Q )-A 



If d 9l T^(A, 9 U . . . , 9 k ) = {uh g ) 2 (9~)/2 - {u\ fi ) 2 {9J)/2 = then u L A e is a true solution of the quasi- 
periodically forced pendulum (3.12). As in the previous section the variation in the action variables is 
given by the partial derivative with respect to A, that is 

d A rf(A, 9) = / + -Jd v f(wt + A, <£(t) + u% e (t)) - d v f(ut + A, q A (t))) dt - (l' Q - I Q ). (3.15) 
J —00 



Lemma 3.4 Let (A, 9) be a critical point of T^. Then there exists a heteroclinic orbit connecting the 
tori Tf and Tjt . 

Id I 

PROOF. By ( 3.15Q it is easy to verify that the solutions of ( |3~l| ) ult + A,q A ' + u L A g , q 1 ^ + u A e ), 

with I„(t) = C-fx f* d v f(us + A, q^s) + u L Ae {s)) ds and C = I' + a„(A) + n / + °° d v f(uit + A, q^(t) + 
u le(*)) - <VM + A, q*X(t)) dt, is a heteroclinic solution connecting 7^ and Tj, . ■ 



Finally, arguing as in the proof of theorem 2.3, we obtain 



Theorem 3.2 Assume (HI) and let satisfy the "splitting condition" 2.1. Then \/Io,I' with u> ■ Iq — 
u> ■ I'q, there is a heteroclinic orbit connecting the invariant tori Tj 1 and Tp . The same estimate on the 

holds. 



'usion time given in theorem 2. 



A Taylor expansion in fi gives 



Lemma 3.5 For \i small enough 

Gp(A) = const + fj,M(A) + 0(i-i 2 ), VAeT" (3.16) 



where M(A) is the Poincare-Melnikov primitive M(A) = f{wt + A, qo(t)) — f{ojt + A, 0) 



dt. 



Proof. We develop with a Taylor expansion in /i the Lagrangian defined in ( 3.13 ) 



_ u 2 



+ (1 - cos u) + n ((u - sin u) 7 + f(wt+ A, u) - f(wt + A, 0) - d g f(wt + A, 0)uj + K(n, u, t) (3.17) 



where 7 (t) := dp, =0 <Z^(*); I^(m>' u )*)I = (^t 2 ), 7?.(/^, 0, i) = and d u TZ(fi, 0, t) = 0. The Mclnikov function 



corresponding to Lagrangian ( 3.17 ) is 



M*(A)= / (q o (t)-siaq o (t))i(t) + f(ut + A,q o (t))-f(ut + A,0)-dJ(ut + A,0)q o (t)dt. (3.18) 

JR 

Integrating by parts, since —7 + 7 = d q f(uit + A, 0), we have 

(qo(t)-sinq (t))ry(t)dt= [ (q (t)-q (t))^(t) dt = [ (-^(t)+^(t))q (t)dt = [ d q f(ut+A,0)q (t)dt, 
r Jr Jr Jtl 



and we deduce from ( gig ) that M*(A) = M(A) = J R [/M + A, q (t)) - /(wi + A, 0)] dt. 



Theorem 3.3 Assume (HI) and let M possess a proper minimum (or maximum) Ag, i .e. 3r > such 
that m{QB r (A ) r > r(^4o)- Then, for /i small enough, the same statement of theorem Iff. 4 feoZds where the 
fusion time is Td = 0((1/ ji) log(l//i)). 



Remark 3.1 By theorems 3.1-S.i. we obtain that, for a priori-stable, isochronous, degenerate systems 
considered in \lu^ 

„2 

U e = euj ■ I + y + e d (cosq - 1) + fif(ip, q) with 1 < d < 2, 

for [i — Se d , 5 being a small constant, the diffusion time is bounded by Td — 0(C(S)/e d ). This improves 
the result of which holds for [i = 0(e d ) , d! > d/2 + 3, and provides the upper bound on the diffusion 
time T d = 0(l/e c+2( - T+1) ( 2d '~ 1 ~ d / 2) ), C is a suitable positive constant. 

4 Splitting of separatrices 

If the frequency vector lo = ui £ contains some "fast frequencies" = 0(l/e b ), b > 0, e being a small 
parameter, and if the perturbation is analytical, the oscillations of the Melnikov function along some 



directions turn out to be exponentially small with respect to e. Hence the development (3.16) will 
provide a valid measure of the splitting only for fi exponentially small with respect to e. In order to 
justify the dominance of the Poincare-Melnikov function when fi = 0(e p ) we need more refined estimates 
for the error. The classical way to overcome this difficulty would be to extend analytically the function 
Ffj,(A, 6) for complex values of the variables, see j§|-|13 and psfl . However it turns out that the function 



F fl (A, 9) can not be easily analytically extended in a sufficiently wide complex strip (roughly speaking, the 
condition q^ g (Re 6) = ir appearing naturally when we try to extend the definition of q^ g to £ C breaks 
analyticity) . We bypass this problem considering the action functional evaluated on different "1-bump 
pseudo-homoclinic solutions" Q^g- This new "reduced action functional" F^(A,0) = G ^{A + wff) has 
the advantage to have an analytical extension in (A, 9) in a wide complex strip. Moreover we will show 



that the homoclinic functions G M , G M corresponding to both reductions are the same up to a change of 
variables of the torus close to the identity. This enables to recover enough information on the homoclinic 
function G^ to construct diffusion orbits. 

We assume that f{f,q) = (1 — cosq)f(<p), f(cp) = J2kez™ /fc ex P(*^ ' f) an d that, there are ri > 
such that 

Vs € N, 3C S > such that |/ fc | < — L exp ( - r *\ k *\) > Vfc G Z ™- t 4 - 1 ) 



Condition (4.1) means that / has a G°° extension defined in 

D := (R + i[-n,n]) x ... x (R + i[-r n ,r n }) 

which is holomorphic w.r.t. the variables for which r, > in (R + x . . . x (R + il n ), where Ii = {0} 
if Ti = 0, Ii = (— fj,ri) if > 0. We denote the supremum of |/| over D as 

:= sup \f(A)\. (4.2) 



It will be used from subsection 4.2 



4.1 The change of coordinates 

Define ip ■ R -> R by tp (t) = cosh 2 (t)/(l + cosht) 3 and set ijjg(t) = tp(t-9). Note that J R Vo(*)9o(*) d< 



7^0. Arguing as in lemma 2.1 



we can prove 



Lemma 4.1 For fi small enough (independently of u>), W € R, there exists a unique function Q A g (t) 
R — > R, and a constant a A g smooth in (A, 9, /z), swc/i i/iai 

• W-Ql 9 W+sinQ^(t)=/i BinQ%{t)f{Lot + A) + a^gMt); 

• («) IniQlAt) - le(t))Mt) dt = 0; 

. {iii)m^Q% e {t)-qe{t)\^Q%{t)-q e {t)\) = o(/zexp(-^l)); 
Moreover Q Afj {t) = Q^ A+k2 ^ e (t), Vfc € Z" and Q^ fl+ ,(t + 7?) = Q^+^W, V0, ?7 G R. 



We define the function F^(A, 9) : T™ xR^Ras the action function al o f Lagrangian 2J3 evaluated 
on the "1-bump pseudo-homoclinic solutions" Q A g (t) obtained in lemma 11, namely 



F^(A,9):= [ C^Q% e {t),Q% 6 {t),t) dt (4.3) 
Jr 

and G^(A) : T" -> R as G M (A) = 0). The following invariance property holds + rf) = 

Ffj, ( A + cur], 9) , V0, jj € R; in particular F M (A, 0) = G M (A + w0) , V0 € R. 



Remark 4.1 By lemma (.l-(i)-(ii), if dgF^(A,9) = then Q Ag is a true solution of (2A). More 

precisely we have \a. Ag \ < C \dgF u (A, 9)\, for a suitable positive contant G > 0. In addition we could 
easily prove using lemma that 

\V'G^A)\ = 0(/x), \V S F»(A,9)\ =0(ji), 8 = 1,2. (4.4) 

The relation between the two functions F^(A,6) and F^(A,6) is given below. The next theo rem is 
formulated to handle with also non-analytical perturbations /. For the analytical case see remark ["" 



Theorem 4.1 For [i small enough (independently of u) there exists a Lipschitz continuous function 
VT n xR^R, wtihh ll (A,6)=0(ji), \h»(A\0j-h»(A,6)\ = 0(^(\A , -A\+\6 / -d\)),h f (A,6+r)) = 
hfj,(A + rju>,6), such that F^(A,9) = F^{A,9 + h^(A,9)). In particular, setting g^A) — h^{A,0), the 
homeomorphism ip^ : T" — > T™ given by t/j^ (A) = A + g^ (A)ui satisfies = G M o tp^ . 



In order to prove theorem LI wc need the next two lemmas, proved in the appendix. 

Lemma 4.2 For fi small enough (independently of u>) there exists a smooth function 1^{A,9) with 
l^A, 9) = 0(p), Vl^{A, 9) = 00), 9 + v ) = l^A + r?o;, 9) such that Q%{9 + l^A, 9)) = n. 

Define V„(A, 9) = F^A, 9 + l^A, 9)). 

Lemma 4.3 There exists a positive constant C4 such that, for all (A, 9) G T™ x R, there holds 

I F M ( A, 9) - V„ ( A, 9) I < C 4 1 do F» ( A, 9) \ 2 . 

In particular if dgF^A, 6)=0 then F^A, 9) = V^A, 9). 



PROOF of Theorem |4.l| . By lemma |12|, there is a smooth function l^ such that l^A, 9) = O(^), 
Vlf,(A,9) = 0{n) l^{A,e + ri) = l^(A + r)uj,9) and F^(A,9) = V^(A, 6 + l^(A, 6)). So it is enough to find 
h = h^(A, 9) such that 

V li {A,e) = F tl {A,6 + h), (4.5) 

h fi (A, 9) will be then defined by 

h^{A,6)= 1^,6) + h^(A,e + 1^,6)). (4.6) 



Note that if deF ll (A, 9) = then, by lemma 4.3, equation (4.5) is solved by h = 0. In general we look for 
h of the form h = dgF^(A, 9)g. Then we can write 

F^(A,9 + h) = F^A,9)+d e F^A,9)h + R^A,9,h)h 2 (4.7) 
= F^A, 9) + (dgF^A, 9)) 2 g + R^A, 9, dgF^A, 6)g){d e %{A, 9)) 2 g 2 

where 

R^A, 0,h) = ^ [F^A, e + h)- F^A, 9) - dgF^A, 9)h 



is smooth and, by t he estimates^ 4.4) on the derivatives of F^, it satisfies R^{A, 9, h) — 0(/i), dhR^{A, 9, h) 
0(/x/|/i|). By ( |4.7|) equation (|4.5|) is then equivalent to 



{d e F^A,9)f 



g + R ti (A 1 9,d e F tl (A,9)g)g 2 



We have R^A,9,d e F^A,9)g)g 2 = 0(ng 2 ) and d g ^{A,9,d g F f2 (A,9)g)g 2 ^ = 0{^g). By the contrac- 
tion mapping theorem, for /i small enough, for all y £ R such that \y\ < 2C4, there exists a unique 
solution g — <p(fj,, A, 9, y) of the equation 

y = g + R^A, 9, d e F^(A, 9)g)g 2 , (4.8) 

such that |<?| < 3C4. Moreover, the function ip defined in this way is smooth. Setting 

h,(A, 6) := <pU A, 9, VM fi~ F it ,9) ) d ^ ^ 
v (dgF^A, 9)) 2 ' 

if dgF^{A,9) 7^ and h^(A,9) = if deF^(A,9) — 0, we get a continuous function h^ which satisfies 
(fy| and \h^(A,9)\ < 3C 4 |<9 e F M (A, 9)\, which implies = 0{n). Moreover h M is the unique function 
that enjoys this properties. 



By ( |4.9| ) the restriction of to 

U„ := {(A, 8) e T" x R : ^F M (i4, 0) ^ 0} 

is smooth. Deriving the identity V fl (A, 8) — F^(A, 8 + h fJi (A, 8)) we obtain 

deF^A, 8 + h^{A, 8))Vh^(A, 8)) = VV^A, 8) - VF^{A, 8 + h^A, 8)) 

for (A,8) e Up. Now, from 

\F„(A,6)-V„(A,6)\ <C 4 \d e F^A,8)\ 2 , d 2 MA,9) = 0(jt), d 2 ee F„(A, 6) = O(n), 

we can derive by an elementary argument whose main ingredient is Taylor formula that 

\VF fl (A > 8)-VVp(A > 8)\ =0(^Ji\deF^A,e)\). (4.10) 

Mo reove r, by the estimates of |<9^F M | and {d g F^){A, 8 + hp(A,6)) = (1 + O^deF^A.O). Hence 
by ( 4 . 1 0| ) |V/i M | = 0{y/~jl) uniformly in Since /i M is continuous and h^A^d) — if (A, 8) £ U^, the 



Lipschitz continuity of follows. 

To complete the pr oof, we observe that hfj,(A,8 + rf) = h^(A + nui,8), which is a consequence of 
uniqueness. Hence, by ( |4.6| ) and the properties of 1^, satisfies the same. ■ 



Remark 4.2 Assume that satisfies the "splitting condition" 2.1 (or its generalisation introduced in 
remark 2.1), with bounds 8 and a. Then, by theorem too satisfies this condition, with constants 



8' , a' . We can take 8' = 8; moreover, at least if -»/jtZ|cj| is small, we can take a' = a/2. As a consequence, 
the results that we shall obtain in the next section proving a "splitting condition" for G^, may be used to 



apply the shadowing theorem 2.3 



Remark 4.3 Assume that > for all i (i.e. that the perturbation f is analytical). Then we can 
prove, using the arguments of the next subsection, that the homoclinic function G M (-) = -Fjj(-,0) can be 
extended to a complex analytical function over the interior of D. Hence F^(A,8) — G^(A + uj9) can be 
defined in an open neighbourhood o/T™ X R in (T n + iR") x C, so that the extension is analytical. One 



could check that l^ and V^, defined in lemma J^.i have analytical extensions too, and that the inequality 



of lemma \4.q still holds in the new set of definition. Moreover in the next lemma 4-4 it is proved that 
F^ is analytical w.r.t. (A,8). As a consequence, (V^A^O) — F^(A, 9))/(dgF fi (A, 8)) 2 is real analytical, 
and so is the function defined in the proof of theorem 4-1- Therefore if Vi > for all i, then the 



homeomorphism tpn defined in theorem 4-1 is a real analytical diffeomorphism. 



4.2 Analytical extension 

The unpertubed homoclinic qa(t) = 4arctg e_ ca n be extended to a holomorphic function over the strip 
S := R + i{— 7r/2, 7r/2). More over equation (2.2) may be considered also for complex values of q and, for 
/i = 0, qg is a solution of (2^2) for all 8 € S. We shall use the notation S$ — R + i(— (f — b), \ — 8) , for 
8 G (0,7r/2). We have 

2 sinh(,z — 9) 2 
qe{z) = —m — 7TT> 1s{z) = smqe(z) = -2 -j- — , (1 - cos qe[z)) = -j- 



cosh(z-fl)' " wvy cosh 2 (z-0)' v ~~*»v~" cosh 2 (z-#)' 

Assume that 8 £ Ss, Re{8) — 0. The following estimates hold, where ieR. 

^ . TtS-IR n exp(-|t|); (4.11) 
min{(|i| + b), 1} 

|sinfffl(*)l < • S (u\,xv 1X exp(-|t|); (4.12) 
min{(|i| + by, 1} 

|C ° S9eWI " imn{(|t|+^,l} ; (4 ' 13) 

' < Cexp(|t|)min{(|t| +8),1}. (4.14) 



Mt)\ 



In what follows we consider the Banach spaces 

! = C 2 (R,C) supexp(|t|/2)(Kt)| + |i&(t)| + \w(t)\) < +00} 



and 



endowed with norm 



X 



{weX w(0) = o\, 



\\w\\ 2 , s = sup (\w(t)\ + \w(t)\ + \w(t)\) exp(H) + sup (-^L + W^L + H t)| 

i*i>i v ' 2 i*i<i Mr + s r (1*1+0) 



\w(t)\ 



The next lemma extends lemma 4.1 for complex values of the variables. First note that the function 
ipo(t) can be extended to a holomorphic function on R + i(— tt, 7f). Recalling the definition for ||/|| given 
in (f4.2|), we have 



Lemma 4.4 There exist positive constants 77, C5 such that for all 5 G (0, tt/2), V0 < [i < (?/^ 3 )/||/||, for 
all to , for all A S D, for all 9 € Ss there exist a unique Q A g : R — > C anrf a unique a A 9 £ C swc/i that 

• #1,0 = 90+^, + w %> where v A,e & C, w € X and \\w%\\ 2 ,s + W%\ + \ a %\ ^ C 5^\\f\\/S 2 ; 

• / R (Qii, 9 (*)-?fl(*))^(*) * = o. 

Moreover Q A g and a A g depend analytically on 9 and on the Ai for which > 0. 
Proof. 1st step. Let us consider the Banach space 

1*1 ■ 



r = {weC(R,C) sup|v(i)|exp(LT) < +00} 



endowed with norm = sup| t | >:l \v(t)\ exp(Jj) +supui <1 (|t| + Let 9 € Ss be given once for 

all. We may assume without loss of generality that Re(9) = 0. 

For 9' S Sg/2 such that \9' — 0\ < 8/2 we introduce the linear operator Lgi : X x C —tY defined by 

Lgi (w,a) = — w + (cos qe>)w — aipg . 

Using that qgi is a solution of — ij + cosqgiy — we can compute the inverse of Lg*. It is given by 
L g , (g) — (w, a) with 

a = _ S».9(t)4r(t)dt 
f R i>e (t)qe'(t) dt 



w(t) 



qe'(t) 
qe>{t) 



( / (.9( cr ) +ail'g(a))qg l (a) daj 



ds 



%>{s) 



+00 



(g(a) + atpg{a))qgi(a) da) ds 



(4.16) 
(4.17) 



Note that since \6-9'\ < 5/2, Re(9') < 6/2. Therefore estimat es ([4 .ll|)-(|4 14[) hold as well (with perhaps 
different constants) when 9 is replaced by 9' . Wc derive from ( 4.15 )-( 4.17 ) that 



M + IM|a, 4 < j\\9\\-i,s. 



(4.18) 



2nd Step. We shall search Q as Q = qe+ v + w with \v\ < 5/2, w S X. Let B denote the open ball 
of radius 5/2 in C centered at 0. Let : B x X x C —> Y x C be defined by 

J^v, w, a) = - ?D + sin(g e+l/ + w) - [isin(q e+u + w) f(tp) - aipg, / (q g+u + w- q e )tpe(t)j- 



From now we shall use the norms \\(i/, w, ct) 1 1 2 = \v\ + IMh.i + \a\ on B X X x C and 
||sil-i,<5 + on y x C. J M is of class C 1 and 

DJ^{v,w,a)[z,W,a] = ( - z + cos(g e+l/ + w)(zq e+u + W) 

- ficos(q e+l , + w){zq e+v + W) /(</?) - a^e, / + W)ip e (t) 



R 



We shall prove that, provided w, a)\\2/5 and yu||/||/i5 3 are small enough DJ^v^w^a) is invertible. 
We first consider the case when w = and \i = 0. Let T„ = DJ (v, 0, a) (independent of a). Observing 

that — <7 e+ „ + cos{qg +v )qe +l/ = 0, we obtain 

T v [z,W,a] = ( - W + cos q e+v W - aipe, J {zqe+u + W)ipe(t)) . 
Using the first step we derive that T v is invertible and, for a suitable positive constant C 

\\T-\g,m2<^\\{9,m-i- (4-19) 
Now we estimate \\{DJ^(i>, w, a) — T v )[z, W, a]||-i- We have 

(DJ^WjO) - T v )[z, W,a] = ((cos(q e+ly + w) - (cosq e+v )){zq e+l , + W) 

- (J,cos(q 0+U + w)(zq e+l , + W)f(tp),0 

We easily get 

\\(DJ^w,a)-T v )[z,W,a]\U < C\\w\\ 2 , g (\\W\\ 2 ,6 + \z\) + ^\z\ + |wn./ linn,, 

< c(|M| 2itf + ^li)||(z,w;a)|| 3 . 



IKDJ^w,^)-^ < 



As a consequence, recalling ( 4.19 ), if //||/||/(5 3 < K and ||io||2,a/<5 < Kq, for if small enough, then 
DJuQi, w, a) is invertible and 

Ki 
5 

for a suitable positive constant K\. 

3rd Step. We now prove the existence of a constant K 2 such that (0, 0, 0) is the unique solution of 
the equation Jq(v,w, a) = in B{K25), ball centered at the origin and of radius K25 for the norm || 
First we observe that, since qe+v = sin(gg +i ,), there holds 



J Q (i/,w,a) = T v \v, w, a] + (sin(g e+ „ + w) - sin(q e+l ,) - cos(q e+u )w, / (q 9+v - qe - fqe+u)^e J ■ 

<j R 

Moreover, by the analyticity of qa,qo, "00 over S, there holds 

(Qe+u(t) - qe(t) - vqe+u(t))ipe{t) dt= (q u (t) - q (t) - vq u {t))%(t) dt, 
r Jr 

hence there is a constant C such that 

(sm.{q e+u +w)- sm(q e+u ) - cos{q g+u )w, J (q e+v - q e - vqe+ v )^e) i < ^'(IIHIl.a + W\ 2 )- 



So, if Jq{v, to, a) = then, by ( |4.19| ) 

- T" 1 ( sin(ge +1/ + to) - sin(g e+ „) - cos(g0+„)to, / (g e+1/ - qg - uq g+l/ )ipe 



(v, w,a) 



< 



CC 



Let K 2 < 1/(CC), By the latter inequality, if Jq(v, w,a) =0 and \ to, ck)| |a < K28, then v = 0, to = 0, 
a = 0. 



4th step. By the previous steps we know that there exist positive constants Kq, Ki and K2 such 



that 



(i) ( Jo(f, w,a) = and | w, a)\ | 2 < K 2 S) 



0: 



(ii) lf\v\<S/2,\\w\\2,s<K S,fi 
Kt/6. 



< K S 3 then DJ^z/, to, a) is invert ible and ||(DJ At (^, u>,a)) 1 || < 



Moreover there exists a constant if a > such that 



(Hi) I {d^J^v, w, a) 



(ain(q e+u + w)f(ip),0) 



< 



K 3 /6. 



We say that (i),(ii),(iii) imply that there is ij such that, for all < fi < ?/<5 3 /||/||, the equation 
J^ L (v,w,a) — has a unique solution such that \\(v, w, ck) 1 12 < K26/2. In addition \\(y, to, ot)\\2 — 
0(^\\f\\/8 2 ). To prove existence, we can proceed as follows. Let S denote the set of all \i G [0, ATo<5 3 /| |/| |] 
such that there exists a C 1 function 5^ : [0, /j] — > {(v, w, a) : \\(v,w, a)|| 2 < K 2 5/2} such that S^(0) = 0, 
Jt{S^,(t)) — for all t e [0,(J>]. S is a bounded interval. Let us call ~p its supremum. By (ii) and the 
Implicit Function Theorem, ~p, > 0. In addition, for fi £ S, there is a unique function with the 
required properties. As a consequence, for < fj, < fi' , 5 M = *S^ M ' | ro,/x] an( i we can define a C 1 function 
5 : [0,~p) {(i/,to,a) : ||(^, to, a)|| 2 < K 2 S/2} such that = S^t) for all ^ € (0,/Z). By (ii) and 
(m), we can write, for all t G (0,/H), 



Hence 



||S'(t)|| 2 = [DJ t (S(t)) 

\\S(t)\\2< 



1 

A 

<5 2 



< 



K X K Z 



.20) 



Now, since S'(t) is bounded, S(t) converges to some S as t — > /I. Either /I = -Ko<5 3 /||/|| or ||5||2 = K28/2 
(If not, by the Implicit Function Theorem, we could exte nd th e solution S to an interval [0, /!+£), £ > 0, 
contradicting the definition of /!)• In the latter case, by ( 4.20 ), 



K 2 8 K 1 K S _ 



S 



So the existence assertion holds for < fj, < ?]S 3 /\\f\\, where i] = min(K , K 2 /(2KiK 3 j). 

In order to prove uniqueness, we assume that there are 61,62 such that ||6i||2 < K 2 S/2, J^(bi) = 0. 
Then, by the same argument as previously, we can prove the existence of two functions of class C 1 
Si,S 2 : [0,/x] -* {b ■ \\b\\ 2 < K 2 S} such that S t (fi) = b u J t (Si(t)) = 0. Moreover, by (ii) and the 
Implicit Function Theorem, Si(^) ^ S2 (/i) implies that Si(t) ^ S 2 (i) for all t € [0, /x], which contradicts 
(i), proving uniqueness. 



The bound of 



u A,e 



A. 9 1 



e I given in the statement is a direct consequence of ( |L20| ). 



To complete the proof, we point out that J p is analytical on (A, 9). Therefore, as a consequence of 
the Implicit Function Theorem (see for example [0), 



Ta,b — 



9+u*{Afi) +W^(A,( 



depends analytically on 9 and on Ai if > 0. ■ 



We now consider the analytical extension of the function F^(A, 8) for (A, 9) G D x Sg 
~ f (& ) 2 (t) 

%{A, 0)= / A ' 6J + (1 - cosQ^(t)) + fx(co S Q%(t) - l)f(ut + A) dt. 

•J R. 

Let consider also the analytical extension for (A, 9) E D x 5,5 of the Melnikov function 

M(A,6) = / (1 - cosq e (t)))f(ujt + A) dt. 
Jr 

We have T(A + uj9) = M (A, 9). We now prove 

Lemma 4.5 For /i||/||i5~ 3 small enough, for all (A, 8) G D x Ss, we have 

,211 f\\2. 



F^A, 9) = const + (j,M(A, 9) + o(^^-) . (4.21) 



Proof. We have g — q g +„^ + g and we set for brevity Q\ g = q g+u + w. 

F„{A, 6) = / qe+v W + {1 - cos{q e+lJ + w))+fi(cos{q e+u + w)~l)f{ujt + A) dt 

= const + I -q g+v w + ^-w 2 + (cosq g+ly - cos(q g+ly + w)) 
Jr 1 

+ pt(l — cos(q g+v )f{ujt + A) + ^(cos q g+lJ - cos(q g+v + w))f(u)t + A) 

= const + fj,M(6 + v, A) + J —w 2 + ^ cos q g + v — cos(q g+l/ + w) — sin q g + v w^j 

+ (i J (^cosq g+u ~ cos(ge+iy + w)jf(ut + A). 

By the estimate ||ia>||2,<s _! CVI l/l I/<5 2 j it follows easily 

Fp (9,A)= Const + »M(6 + v, A) + O ( ^ |j| 1 1 - ) . 

For example we can get that J R cosge+i/ — cos(<7e+iy + — (shige+i/)^ = 0(/i 2 ||/|| 2 /<5 ) by writing 
cos q g+v — cos{q g+u + w) — (sin q g+u )w — w 2 f Q — (1 — s) cos(q g+iy + sw) ds and using ( |4. 12| )- ( 4. 15 ) togheter 
with | |if 1 1 2,5 < /x||/||/5 2 . Moreover 



\M(9 + v,A)- M(9, A)\=0(^)= , 



which completes the proof of the lemma. ■ 

The Fourier coefficients of the Melnikov function T(A) = J^k ^ k ex p(*fc^4) are explicitely given by 

2n(k-uj) 
smh((/e • w)§) 



By estimate ( |h2l| ), since F^A, 8) = G^(A + u9) and M(A, 9) = T(A + uj9), via a standard lemma on 
Fourier coefficients of analytical functions (lemma 3 in we obtain the following result (compare with 
theorem 3.4.5 in §(J). 

Theorem 4.2 There exists a positive constant Cq such that, for /x||/||<5 -3 small enough, then\fk ^ 0, k € 
Z n , for all S £ (0, §), for alluj, 



G fc -Mr fc | < ^1' " exp ( - J2 r i\ki\) exp ( - \k ■ u\ (| - s)) ■ (4.23) 



5 Three time scales 

We consider in this section three time scales systems as (see jy| and [p2| ) 

h p 2 
U = -±= + e a (3 ■ I 2 + — + (cosg - 1) + //(cosg - l)f((p. u <p 2 ), e>0 

with n > 2, ipi £ T 1 , tp 2 £ T™ -1 , I\ £ R 1 , 7 2 £ R n_1 , (3 £ R™ -1 and e is a positive small parameter. 
The frequency vector is uj = (1/^/e, e a /3), where (3 — (/3 2 , ■ ■ ■ , Pn) G R™ -1 is given. 
We assume through this section that p\ |/| |e~ 3 / 2 and e are small. 

Given k 2 = (fe, ■ ■ ■ , k n ) £ Z n_1 , we shall use the notation n\ := (|A^j|, . . . , \k n \). Moreover we shall use 
the abbrevi atio n p 2 := (r 2 , . . . , r n ), so that n\ ■ p2 := X]"=2 r »IM- We recall that r\, . . . , r n are defined 
in formula ( |4.l| ). 
Writing 

f{<Pi,<P2)= ^2 fk 1 .K 2 exp(i(k 1 ip 1 + K 2 -(p2)), 
(fei,K 2 )eZxZ"- 1 

we assume that / is analytical w.r.t tf2- More precisely, 7*1 = and for i > 2, > 0. If a = 0, we impose 
in addition that > |/3j|7r/2 for i > 2. 



We shall use (4.23) in order to give an expansion for the "homoclinic function" 



Gp{A) = Gk ^ exp(i(fc 1 A 1 + k 2 • A 2 )) = ^ 9k x {M) exp(ifci^i). 

(fci,K2)eZxZ"- 1 fciez 

We start with 

Lemma 5.1 There exists a positive constant C7 such that, for /z||/||£ -3 / 2 small enough, 



E 



K 2 eZ"- 1 ,|fei|>2 



|G fcl , K2 |<C 7 ^exp(-^). 



(5.1) 



PROOF. Choosing 8 = ^fe, we get from (4.22) and (4.23) 



ki 



< 



Cp\\f\\e-<^( 



k 2 ■ (3s a 



I \ e -\^+^2-f}^\n/2 



+ C^_||J||2 e -Kj-P2 e -|^ + «2-/3£ a |(7r/2-VF) 

< C^Pdfcil + | K2 |) e -4-p 2 +l«=^l^-/2 e -^-/2 

+ (7^||f||2 e -K+-P2 + |K2^k a ( 7 r/2-v^) e -|fcll( 3 7 ; -l) 



< 



cMml - n(\h\ + |« 2 |)exp(- £ |^|(r, - 10^/2)) exp(-|fc 1 |(-^ - 1)). 

V £ j=2 ly J £ 



We have used in the last line that /x||/||/e 3 / 2 = 0(1). Now r, — \(3j\e a n/2 > for e small enough both if 
a = or if a > 0. Summing in |fci| > 2 and in k 2 € Z™ -1 we obtain ( |5.l| ). ■ 

Let 



r(e,A)= ^ r fcliK2 exp(i(fciA 1 + k 2 • A 2 )) = ^ r fel (£,A 2 )exp(i/ciAi). 
(fci.^leZxZ"- 1 feiGZ 



Lemma 5.2 We have 



g {M) = Me,A 2 ) + O(i?\\f\\ 2 ). 



Proof. A summation over k 2 in estimate ( 4.23 ) (where we chose S — ir/2 and k\ = 0) yields immediately 
the estimate. ■ 



Lemma 5.3 We have 



9±i(A 2 ) = fxT ±1 (e,A 2 ) + (O^L exp(-^=)) 



PROOF. By ( 4.23| ) (where we chose 5 — \fe and k\ — ±1), we can obtain as in the proof of lemma 5J 



9±i(A 2 )-vT ±1 (e,A 2 )\ < C^||/|| 2 J2 exp(-^fe|(r J --| ( 5 i | £ V2))exp(-(^=-l). 



< ^ 



Since T(A) and G M (A) are real functions we have that g-i(A 2 ) = g 1 (A 2 ) and T_i(A 2 ) = Ti(A 2 ), 
where z denotes the complex conjugate of the complex number z. We deduce from the previous three 
lemmas the following result. 

Theorem 5.1 For (J-\\f\\s~ 3 / 2 small there holds 

G^(A 1 ,A 2 ) = Const + ^fiT (e,A2)+Ro(e,n,A 2 )j + 2Re ^iT^e, A 2 ) + R^e, a, A 2 )^e lAl 



0(M£- 1/2 ||/|| exp(-— )) 



whe 



R (e, \i,A 2 ) = 0[n 



and R\(e, \i, A 2 ) = O 



■ exp(- 



Remark 5.1 (i) This improves the results in 111/ which require /i — e p with p > 2 + a. 



(ii) The 



5.1 certainly holds in any dimension, while the results of Uaj, which hold for more 



general systems, are proved for 2 rotators only. 



(Hi) Theorem 5.1 is not in contradiction with the counterexample given in 

(iv) In order to prove a splitting condition using theorem 5.1 it is necessary, according with [l]\ ] and 
J2lfy, that 3m, I G Z" _1 such that /o,;,/i,m =/= 0. Otherwise, recalling ((.2i), it results that Tq(e,A 2 ) 
E« 2 ez«-i rV 2 exp 4K2 - A * = and also Ti(s, A 2 ) = £ K26 z«-i r 1)Ha exp*" = 0. 



Theorem 5A enables us to provide conditions implying the existence of diffusion orbits. For instance 
we obtain the following result. 

Lemma 5.4 Assume that there are A 2 € R" _1 and do,co > such that, for all small e > 0, 

(i) |ri(e, A 2 )\ > {co/y/e)e-*t<*ft, VA 2 e R"" 1 such that \A 2 _-A 2 \ < d ; 
(ii) T (e,A 2 ) > T (e,A 2 ) + c , VA 2 e R' i_1 , such that \A 2 -A 2 \=d . 



Then there is C\ > such that, for ju||/||e 3 / 2 small enough, condition 2.1 is satisfied by G„, with 
a = ^e- 7 ^ 2 ^) /y/i and S = Co/V^x/i^" 7 ^ 2 ^. 



Proof. First we can derive from ( 4.22| ) and (4.23) in the same way as in the proof of lemmas 5.1 and 
5.3 that 

|ft(4,)| + |Vft(i4 2 )| < Yl (i + M)|G ljK2 | 

K2GZ"- 1 

(5.2) 



< c M\J\\ c -n/(2V?) 



By the bounds of R and R± of theorem 5.1, for e and /i||/||e 3 ^ 2 small enough, we have 

(i) \gi(A*)\ = K/xTi + Ri)(A 2 )\ > ( MC o/(2^)_)_e-" /(2 ^ ) VA 2 e B do 
(ii) |?o(A 2 )| = (^r + i?o)(^ 2 ) > (^r + Ra)(A 2 ) + co/2 VA 2 e dB dn , 

where B do is the open ball centered at A 2 of radius do. 

So we can write g\ (A 2 ) — \gi{A 2 )\e % ^ Azi \ where is a smooth function defined in Bd a ( ( |5.2| ) and the 
previous lower bound of |^i (^2) | provide a bound of V<fi(A 2 )). 

For A 2 £ Bd , by theorem 5.1 we have 

G^(Ai,A 2 ) — Const+(fiT Q +Ro)(e, fi, A 2 )+2\(pTi+Ri)(e, fi, A 2 ) \ cos(Ai+(f)(A 2 ))+0(fj,e^ 1 ^ 2 \\f\ \ exp(— -^=)). 



Let 



U = {(A 1 ,A 2 )eHxH n - i : A 2 eB do , \A, + <j>(A 2 ) - tt| < -}. 



If A 2 S «9i3 do then 



G p (Ai, A 2 ) - G m (tt - 0(A 2 ), A 2 ) > ^ + 0(//e" 



ex P (-^)). 



If \A X + 4>{A 2 ) - tt| = tt/2 then 



G, 



,,{A u A 2 )-G^-K-ct>{A 2 ),A 2 ) = 2\^T l +R l )(A 2 )\+0^e- 1 '\f\\eM-^)) 



M c p -*/(2VS) 



0( Ai£ - 1 /2|| / ||exp(--=)). 



Hence, for e and 3 ^ 2 small enough, 



inf G u > inf G u 



c A*e 



2^1 



Using that |VG,j| = 0(/i), we can easily derive that condition |2.1| (not with a ball B p but the bounded 
open set U, according to remark 2.3) is satisfied with 5 = (co/2)[ie~ 7 '/ 2 ^ 1 
some positive constant c\. ■ 



e, a = C\e 



"l^lJe for 



The condition given in the previous lemma is not easily handable. We now want to provide simpler 
conditions, involving properties of the perturbation /. For A = (Ai, A 2 ) € T 1 x T™ -1 , let 

f{A u A 2 ) = ^ f kuK2 e^p(i(hAi + k 2 ■ A 2 )) = ^ A-i {M) exp(ifciAi). 

(fci,K2)GZxZ'*- 1 feiGZ 

If / is analytical in some domain then also fk 1 (') can be analitically extended is the same domain, as 
Ms) = (1/2-k) f(a,s)e- ik ^ da. 

Theorem 5.2 Assume that f satisfies one of the following conditions: 

(i) a > 0, fo(A 2 ) admits a strict local minimum at the point A 2 and fi(A 2 ) =/= 

(ii) a — 0, fo(A 2 ) admits a strict local minimum at the point A 2 and fi(A 2 + i(7r/2)/3) 7^ 0. 



Then, for all small e such that oj £ = (3e a ) satisfies 



_ | k | T , 7- 

for all Iq,I' with tu e ■ Iq = uj t ■ Iq, there is a heteroclinic orbit connecting the invariant tori 7} and Tj^. 
In addition, for all r\ > small enough the "diffusion time" Td needed to go from a r] -neighbourhood of 
T Io to a ^-neighbourhood ofT^ is O(\I - J£|(v^//i)e 7r /( 2 ^) [(j e )~ 1 (Vee 7r ^ 2 ^) T + | ln/i|] + | ln(??)|). 

Proof. It is enough to prove that, if (i) or (ii) is satisfied, then the condition given in lemma holds. 



The statement is then a direct consequence of theorem 2.3 



We first assume that condition (i) is satisfied. In what follows, the notation u = 0(v) means that 
\u\ < C\v\, where C is a universal constant. We have 



T (e,A 2 ) = 



2 ™ 2 • Pe a M 

/0,K2 e 



sinh(7rK 2 • /3e»/2) 
(4 + 0(e 2a | K2 | 2 ))/o,K 2 e^^ 



4/o(A 2 ) + o( Yl z 2a \K2\ 2 e- K i- p2 ) 
4/o(A 2 ) + 0(s 2a ). 



Moreover 
rx(e,A 2 ) 



E 



27r(£- 1 / 2 + k 2 • (3e a ) 
1 smh((ir /2)(e~ 1 / 2 + k 2 - (3s a )) 

4tt 



A,K2 e 



i«2'^2 



n 2 ez 

47T 
4-7T 



C7»-l » 



ZK2'j42 



-(tt/2) £ 



-1/2 



/i(A 2 )+Of e-^'^de-^/ 2 )" 2 -" 6 " - 1| + |K 2 |e a+ ( 1 / 2 ) e l K2 -' 3 l £a ) 



K 2 eZ"-! 3=2 



i!L e -(T/2) £ - 1/2 



/ 1 (A 2 ) + 0(£ a ) , 



provided e is small enough. It is then clear that condition (i) implies that assumption of lemma |5.4| holds. 
We now assume that condition (ii) is satisfied. As previously, we have 



ri ( £ ,A 2 ) = £ ^ e -(^)( S - 1/ W) (1 + 0( | K2 | £ i/2 )/l ^ e 



E 


4-7T 

-i 


K 2 eZ"- 




47T 


(7r/2) E - 






47T 


(ir/2) e - 







n 

/ 1 (A 2 + <(7r/2)^) + o( X! ^pC-EI^I^-^/ 2 )^'!))!^!^ 



/2 



/i(A 2 + i(7r/2)^) + 0(Vi)l. 



J'=2 



We observe also that, if a = 0, then To(£, A 2 ) is independent of s. It follows easily that condition (ii) 
implies that the assumption of lemma 5.4 holds true. ■ 



Remark 5.2 In many examples, condition (i) or condition (ii) is satisfied. However we need that fn(A 2 ) 
and f\(A 2 ) do not vanish everywhere, see remark 5.1-(iv). 



6 Appendix 



In the proof of the following lemmas we will closely follow the arguments developed in the papers |J-|J 
to which we refer for further details. In the sequel the notation u = O(v) (resp. u — o(v)) will mean 
that there is a constant C (resp. a function e(v)) independent of anything except / such that \u\ < C\v\ 
(resp. \u\ < e(v)\v\ and lim II _ + oe(w) = 0). 

Proof of lemma 2.1. We first assume that 8 = and give the existence proof in [0, +oo). We are 
looking for a solution of (2.5) in the form of q = qo + w with w(0) = and lim f ^ +00 w(t) = 0. The 
function w must satisfy the equation 

— w + w = — ^sin(qo + w) — sinqo — iuj + psm(qo + w)f(u>t + A). 



Let 



and 



X = 



ftu(-) e W 1,oo ([0,+oo)) \\w\\! := Bupmax(|ti>(t)|,|w(t)|)exp(J|-) < +00} 
1 ten * > 



X = («;(•) e L°°([0,+oo)) ||w||o := sup|w(i)|exp(^) < +00}. 
^ teR 2 > 

X and X', endowed respectively with norms || ||i and || ||o, are Banach spaces. Let Cq be the linear 
operator which assigns to h £ X' the unique solution u = CqH of the problem: 

— u + u = h 

u(0) = , lim t ^ +0O u{t) = 0. 



An explicit computation shows that, for t € [0, +00), 

1 f + °° / 

u{t) = (£ h)(t) = -j (t 



-\ts\ _ 



e- ( - t+s Ah(s)ds. 



(6.1) 



As an easy consequence Cq sends X' into X continuously. 

We define the non-linear operator H : R x R™ x X — > X by 



H(p, A, w) := w — £ ( - ( sm(q + w) — sinq - u?J + psin(q + w)f{ujt + A)j . (6.2) 

H is smooth, 27rZ"-periodic w.r.t. A and we have H (0,^4,0) = 0. The unknown w must solve the 
equation H(p,A,w) = 0. We can apply the Implicit Function Theorem. In fact, let us check that 

d w H(0,A,0): W ->W-CoUl-co8q )W 

is invertible. Since lim t _ +00 (l — cos qo(t)) = 0, d w H(0, A, 0) is of the type "Identity + Compact" and then 
it is sufficient to show that it is injective. W is in the kernel of d w H(0, A, 0) iff W(0) = and W satisfies 
in (0, +00) the equation 

- W + cos q W = 0. (6.3) 



Multiplying by q in ( |6.3| ) and integrating over [0, +00) by parts twice we obtain that W(0)qo(0) = 0. 
Since g Q (0) ^ we get also W(0) = and then W = 0. Thus the kernel of d w H(0, A, 0) is reduced to 0, 
and this operator is invertible. We derive by the Implicit Function Theorem that there are po > and 
fiQ > such that, for all \p\ < po, for all A £ R™, the equation H (p, A,w) =0 has a unique solution w 1 ^ 
in X such that ||u>^|| < po- 

Note that po and po may be chosen independent of A (and of uj too) because d w H(0, A, 0) is inde- 
pendent of A and u>, 8^11(0, A, 0) is uniformly bounded, and d w H(p, A, w) (resp. d^H(p, A, w)) tend to 
d w H(0,A,0) (resp. ^#(0,^,0)) as (p,w) -> (0,0) uniformly in (A,u). 

Since H is smooth depends smoothly on p and A and w^ +27rk = u/^ by the 27rZ n -periodicity of 
H w.r.t. A. By the properties of d^H mentioned above, ||w^||i = O(p). 



In a similar way we can prove the existence and unicity of u/ A : (— oo, 0] — > R which satisfies 
analogous properties over the interval (— oo, 0]. We can define <f A by q A Q (t) = qo(t) + w A (t) if t > 0, 
Qao = Qo(t) + w 'a(^) if i < 0. This is the unique function for which (i), (ii) (with 9 = 0) and (Hi) hold. 

If 9 7^ 0, we observe that q satisfies (i) iff 

-(T- e q)" + sm(T-eq) = A* sm(T_ e q)f(ut + A + lj9) 
(T_ 6 q)(0)=iT, 

where T_gq(t) = q(t+9). Hence there is a unique q A g which satisfies (i), (ii), defined by q A g = Tgq A+Lug Q , 
i.e. q A g (t) = q A+ujg (f — 0); (Hi) and (iw) clearly hold. The regularity of q\ g w.r.t. A,9,fi is a 
consequence of the regularity of w A and w' A w.r.t. A and /x. (v) follows from 



d A w 



d w H(fi,A,w») d A H(»,A,w*) 



provided we can justify that []d A H(ix, A, = O(fi), \\u ■ d A H(fi, A,w A )\\i = 0(fj,). The second 

bound (uniform in u>) is not so obvious. We just point out that 

uj ■ d A H(n,A,w'X) = -£ (M sin (9o + w%)—f(wt + A ) 
and that we can use the "regularizing" properties of Co ■ ■ 



Proof of lemma 2.4. We give the proof in the interval [#i,02]- We may assume without loss of 
generality that 9\ = since, by the remark at the end of the proof of lemma 2.1, a translation of the 
time by —9\ amounts to adding lo9\ to A. For si mpl icity of notations, we shall write 62 = 0. 

We are looking for a solution q = q^ g + w of ( |2.2| ) over (0, 6) with w(0) = = 0, where Qg e is the 
following smooth "approximate solution" 

f if 4 6(0,0/2-1), 

9o,e(*)= < r|(*)if t €[0/2-1,0/2 + 1] 
I 2tt + if t 6(0/2 +1,6), 

where 

rj(t) = (1 - R(t - e/2))q^ (t) + R(t - 9/2)(q%(t) + 2tt), 

and R : R -> [0, 1] is a C°° function such that R(s) = if s < -1, i?(s) = 1 if s > 1. Let £ 0: <? be the 
linear operator which assigns to h £ L°°([0, 9]) the unique solution u = £o,gh of the problem: 

-u + u = h . s 

u(0) = , u(0) = 0. 1 j 



An explicit computation shows that for t € [0, 9] the solution u of (6.4) is given by 



1 



u(t) = — / h(s) sinh(s) sinh(0 — t) ds + / h(s) sinh(0 — s) sinh(t) ds 

sinh(0) ' ' 







Note that C .g sends L°°([O,0]) into W 1 ' 00 ^, 0]) (W 2 ' 00 ([0,9]) in fact) and that there is a constant K 
independent of 9 such that HAj.e^lli.oo < ^H^llooi where || Hoc denotes the infty norm in [0,9] and 
IWIr^HWU + llWlU. 

We define the smooth non-linear operator ff fl :RxR"x W 1 ' 00 ^, 0]) -> W llO °([0, 9]) by 

i/ e (^, A, w ) :=w- Cqm ( - ( sm{ql e + w) - q* ofi - w) + (J, sm(q^g + w)f(ut + A)j . 

We immediately remark for further purpose that 

\]d 2 ww H s (»,A,w)[W,W}\\ = 0[\\W\\l). (6.5) 



Moreover, by lemma |2.l| -(z) and the definition of g , || — sin^ g + g + ^sin(</Q g )f(ujt + A)j ||oo = 
O(exp(-0/2)) hence 

\\H 6 (p, A0)||i,oo = O(exp(-0/2)). (6.6) 

lo.e +w is a solution of (^) with the appropriate boundary conditions iff H e (p,, A,w) =0 
We shall show that there exist C,L,~p > such that V0 > L, for all \p\ < JI, for all A and u>, 
d w H e (fi,A,0) is invertible and 

| (^ffVAO))" 1 !! <C. (6.7) 
Since d w H (p,, A, 0) is of the type "Id + Compact", it is enough to prove that 

vw s w rl>0 °([o,^]) ||crf(,i, >i,o)w||i,oo > il|w|| ll0 o. 

We shall just sketch the proof of this assertion (see also lemma 2 of H). Arguing by contradic- 
tion, we assume that there are sequences (/x„) — > 0, (9 n ) — > oo, (A„), (w„), (W n ) such that W„ G 

W 1 ' oo ([0,fl n ]),||W n ||l,oo = l, 

||a !i; JI e "(/i„,A„,0)W„||i, oo ^0. (6.8) 



Let £„ S [0, # n ] be such that m n := max tg r Ol i|W„(t)| = W„(£„). By ( |6.8D and the properties of £o,e, 
||Wn||i,oo = 0(m n ). Hence liminf(m„) > 0. Taking a subsequence, we may assume that (£„) is bounded 
or {9 n — £„) is bounded or ((£„) — > oo and [6 n — £„) — > oo). 

In the first ca se, still up to a subsequence W„ — » IF =^ uniformly in compact subsets of [0, oo). 
Taking limits in (|j| we obtain that W(0) = 0, -W" + cos q W = 0, which contradicts W ^ 0. The 
second case can be dealt with similarly. In the third case, up to a subsequence, W n {- + £„) W ^ 
uniformly in compact subsets of R, with |W(i)| < | W(0) | for all t e R. Taking limits in fl6.8|), we obtain 
that — W + W — over R, which contradicts W ^ bounded. 

From now we shall assume that \p\ < po <~p, 9 > L. Let 

i? e (^, A, w) = H® (/i, A, iu) - i? A, 0) - d w H e (p 1 A, 0)w. 
By the previous assertion, 

H e (}i,A,w) = «■ w = -(9 w F e (M )J 4, ! 0))" 1 ^(M !J 4,0)-(a w F e (M,A,0)) _1 i? e (M !J 4,«;) := F^w). 

We just have to show that F^ is a contraction in some ball B(0,p) C M /1,oo ([0, 6 1 ]). For this, we derive 
from (3.5) and ( |6.6| ) in a standard way that, for all ||u>||i,oo ; ||w/||i,<x> < P, \p\ <~P, 9 > L there holds 

lKAMIIi,oo=O(exp(-0/2)+p 2 ) ; - F^ A (w')\\ = 0(p\\w' - w\\). (6.9) 

We can deduce that F® A is a contraction in B(0,p), with p = Cexp(— 0/2), for some constant C, 
provided that 9 > L, L large enough. Applying the Contraction Mapping Theorem we conclude that 
ther e is a unique solution ||u>^(A, 6*)||i j00 < Cexp(— 9/2) of the equation H A {w) = 0. Note that by 
( |6.9|) uniqueness holds in -6(0, po) for some po > independent of 0. The regularity of the solutions in 
(A,9,p) follows like in S. ■ 



PROOF of lemma \i.2\ . Let us consider the function ff :RxT"xRxR^R defined by 

H(n,A,6,l) = Q Ate (6 + l)-Tt 

The unknown Z M (A, 9) can be implicitely defined by the equation H(p, A, 9, 1) — 0. We have H(Q, A, 9, 0) = 
and 

diH(0,A,e,0) =q (0) ^0. 

Hence by the Implicit function theorem, for p small enough (independently of A, 9, u> because diH and 
d^H are continuous uniformly in A, 9, uS), there exists a unique smooth solution l^A, 9) = O(p) of 



H(/i,A,9,l) — 0. Moreover, by the uniform estimates in A and to that we can obtain for 8a Qa q, 
u ■ d A Q%, there holds \Vl^{A,9)\ = O(p). ■ 

Proof of lemma 4.3. The first step is to prove that 



t - 



Ke+UA^W - QleWl l<?W(A, e) (*) " QaM ) < K \d e F,{A, 0)\ exp(-^), Vt G R. 

(6.10) 

We have g^^+^^.e) = T e+i^A,e)q% ja> Ql,e = T e+h<(A,e)QA>,e" where A' = A + w{6 + l^(A,6)) and 
9' = —l^(A, 9). So it is enough to prove the estimate for w 



lA>,6> ~ Qa',0- 



Note that Q A , e ,(0) = q% i0 (0) = n. So Q A , g , — q a , q A , — qo belong to X and satisfy 
Ufa, A', q\, - g ) = ; H(n,A',Q A , g , - q Q ) = a% e ,Co{^e'), 



where X, H and Co are defined in the proof of lemma 2.1. Therefore 

afX, ie >Ca(ipg>) = d w H(fj,,A',q% fi - q )(Q A , i9 , - q% >0 ) + o{\\Q A , <g , -g^ i0 ||i)- 

Moreover HQ4 , gl — go||i + \\qA' — <7o 1 1 1 = O(p). Hence, by the properties of H mentioned in the proof 
of lemma 2.1 (in particular the fact that 8 W H(0, A', 0) is invertible) we obtain, for /i small enough the 
following bound : 

IIQ^-^olli = °(K^I)- 

Since \a% <e ,\ = 0{\8gF^{A, 9)\) we deduce estimate (|6.10D . 
We can now estimate F^A, 9) - V^A, 9). For q G X let 

G%(q)= f £^A>(q,q,t) dt= f Uq{t)f + (1 - «»(?(*))) + /i(cos(g(f)) - l)f(A' + ut) dt. 

By standard arguments, Q A , : X — > R is smooth. Moreover for < fiQ, 

D 2 g%(q)[w,w}= I w 2 +cos(q)w 2 ~ ncos(q)w 2 f{A / + ujt) dt = 0(\\w\\l). 



By the definition of q A , ((z) in lemma 2.1), we easily obtain with an integration by parts that DQ A , (q A , )w 
for all w € X such that w(Q) — 0. Therefore 

Gl>(Q% fi + ™)=g%(f A ', ) + O{\W\l) 

for all w G X such that u>(0) = 0. Hence since (Q A , g , — q A , )(0) = 0, 

F^A, 9) - V^A, 9) = F^A, 9 + l^A, 9) + 9') - F^A, 9 + l^A, 9)) = F^A' , 9') - F^A* , 0) 
= 0%(Ql',o') - Q%{q%, ) = 0(\\Q% )9 , - &,oll?)- 



We obtain by (|6.10|) that 



for some positive constant C4. 



\F^A,9) - V^(A,9)\ < C A (d e F^(A,6) 
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